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摘要 : 48 (Eriophyton wallichii) 是 中 国 青藏 高 原 地 区 高 山 流 石 滩 特 有 植物 。 为 探讨 第 四 纪 冰 期 气候 波动 
对 青藏 高 原 高 山 流 石 滩 植物 居 群 遗传 结构 以 及 谱系 历史 的 影响 ， 检 测 了 绵 参 20 个 居 群 187 个 个 体 的 核糖 
体内 转录 间隔 区 (17S) 的 序列 变异 ， 共 发 现 19 个 单 倍 型 ， 且 结果 表明 大 部 分 居 群 拥有 独自 的 单 倍 型 。 
AMOVA 分 析 表 明 ， 居 群 间 的 变异 占 了 总 变异 量 的 89. 54% ， 居 群 遗传 分 化 指数 很 高 (Csr =0.863 ，Nsw = 
0.957) ， 但 无 明显 的 谱系 结构 。 我 们 推测 这 种 独特 的 谱系 地 理 结构 可 能 主要 是 由 于 高 山 流 石 浴 独特 的 地 
形 和 环境 及 “孤岛 效应 ”形成 的 ， 与 先前 研究 过 的 青藏 高 原 植 物 不 同 ， 绵 参 呈现 出 高 度 分 化 的 遗传 结构 ， 
常常 呈现 出 “一 个 居 群 ， 一 种 单 倍 型 ”的 模式 ， 这 可 能 是 由 于 第 四 纪 冰 期 及 间 冰 期 尽管 纺 参 的 分 布 范围 
有 扩大 和 收缩 ， 但 始终 由 于 特殊 的 高 山 生 境 和 流 石 滩 基 质 ， 导 致 了 该 物种 长 期 以 来 居 和 群 间 很 少 有 基因 交流 
并 一 直 处 于 隔离 状态 ， 造 成 了 现在 这 样 的 谱系 地 理 结 构 。 
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Phylogeographical Study on Eriophyton wallichii ( Labiatae ) 
from Alpine Scree of Qinghai-Tibetan Plateau 


WANG Xiao-Xiong'^, YUE Ji-Pei! " , SUN Hang’, LI Zhi-Min' *™* 
(1 School of Life Science, Yunnan Normal University, Kunming 650092, China; 2 Key Laboratory of Biodiversity and 
Biogeography, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201 , China) 


Abstract; In this study, we examined the nuclear ribosomal DNA internal transcribed spacer (/TS) sequence varia- 
tion of Eriophyton wallichii, a perennial endemic to alpine scree of the Qinghai-Tibetan Plateau. Nineteen haplotypes 
were detected by analyzing 187 individuals from 20 populations of E. wallichii. Most populations hold the unique 
haplotype, which showed different population with different haplotype. An analysis of molecular variance ( 4AMOVA) 
for populations showed that the genetic variation mainly resided among populations (89.54% ) , the level of differen- 
tiation among populations was very high ( G5, 20. 863; N,, 20.957) but didn’t show a significant phylogeographical 
structure. It is different from foregoing other species reported in Qinghai-Tibetan Plateau that the haplotypes of E. 
wallichii showed a distribution pattern that almost one population hold one unique haplotype, hardly share same halo- 
type among populations. We hypothesized the unique genetic structure of E. wallichii might result from allopatric 
fragmentation in the “islands of alpine" due to the extremely isolated alpine scree and rough topography of the re- 
gion, through the Quaternary climate oscillation dramatically caused population expanding or shrinking. 
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近年 来 对 中 国 喜 马 拉 雅 地 区 及 青藏 高 原 的 
植物 类 群 开展 的 谱系 地 理学 研究 已 取得 了 较 大 进 
展 ， 对 了 解 青藏 高 原 隆 升 以 及 第 四 纪 冰 期 对 其 植 
物 多 样 性 的 形成 及 分 布 格局 的 影响 提供 了 重要 的 
证 据 ， 归 纳 起 来 有 2 种 类 型 。 一 是 在 第 四 纪 冰 期 
一 些 植物 类 群 如 祁 连 圆 柏 (Juniperus przewalskii ) 
(Zhang a 2005), 云南 红豆 杉 ( Taxus wallichi- 
ana) (Gao 等 , 2007) 、 青 海 云 杉 ( Picea crassifo- 
lia) (Meng 55, 2007) , KEH iE (Pedicularis 
longiflora) (Yang $F, 2008), 4& cK dis JE A 
( Metagentiana striata) ( Chen 等 , 2008 ) Ji fr 
毒 (Stellera chamajasmae) (Zhang 等 , 2010) ， 从 
气候 寒冷 的 高 海拔 或 高 原 面向 高 原 东部 及 东南 部 
海拔 较 低 的 地 区 转移 ， 在 此 得 以 存留 和 发 展 ， 因 
而 保存 了 较 多 的 单 倍 型 和 较 高 的 遗传 多 样 性 。 在 
间 冰 期 以 及 冰 后 期 气候 回暖 时 ， 植 物 从 高 原 边 缘 
避难 所 向 高 原 面 或 腹地 回迁 ， 因 受 黄 基 者 效应 及 
瓶颈 效应 影响 ， 仅 部 分 单 倍 型 得 以 在 高 原 面 或 高 
原 腹地 上 固定 下 来 ， 形 成 现今 高 原 面 上 的 居 群 单 
倍 型 较 单一 ， 遗 传 多 样 性 也 较 低 ， 而 在 高 原 边 缘 
海拔 较 低 地 区 的 居 群 遗传 多 样 性 较 高 ， 单 倍 型 也 
较为 多 样 的 格局 。 第 二 种 情况 以 露 蕊 乌 头 (Aco- 
nitum gymnandrum) (Wang 等 , 2009a) 、 银 露 梅 
( Potentilla glabra) (Wang 等 , 2009b) 、 西 川 红 景 
天 (Rhodiola alsia) (Gao 等 , 2009)、 西 藏 圆 柏 
( Juniperus tibetica) ( Opgenoorth 等 , 2010) 、 金 露 
梅 (Potentilla fruticosa) (Li 等, 2010) 、 偏 花 报 
春 (Primula secudiflora) ( Wang 等 , 2008) 、 海 仙 
报 春 (Primula poissonii) ( 宋 敏 舒 等 , 2011) 为 
代表 ， 这些 植物 类 群 在 冰期 时 ， 并 非 所 有 居 群 都 
向 高 原 边 缘 低 海拔 地 区 转移 ， 有 部 分 居 群 在 高 原 
台面 局 部 地 区 存留 和 发 展 ,保留 了 较 多 的 古老 单 
倍 型 和 遗传 多 样 性 ， 在 间 冰 期 及 冰 后 期 气候 回暖 
时 ,植物 从 高 原 面 上 的 原生 地 避难 所 向 周边 扩散 
发 展 ， 形 成 现今 高 原 面 上 部 分 地 区 拥有 和 较 高 遗传 
多 样 性 以 及 较 多 的 单 倍 型 格局 。 因 此 ， 与 高 原 复 
杂 的 地 质 历 史 以 及 气候 条 件 对 应 ， 高 原 上 不 同 的 
植物 类 群 也 相应 的 具有 复杂 多 样 的 进化 历史 。 

迄今 为 止 ， 针 对 中 国 喜 马 拉 雅 地 区 和 青藏 高 
原 植物 的 谱系 地 理学 研究 ， 都 是 以 森林 以 及 高 山 
草 多 植 被 带 的 类 群 为 主 开展 的 ， 而 最 极端 环境 的 
高 山 流 石 滩 生态 系统 中 的 植物 类 群 的 谱系 地 理学 




































































研究 ， 到 目前 为 止 还 未 见报 道 。 高 山 流 石 滩 或 称 
高 山 冰 缘 带 是 指 分 布 海拔 一 般 在 (3 900 ~ ) 4 200 
~5 200 m， 与 冰雪 带 相 邻接 的 生态 系统 ， 是 陆 
地 植物 多 样 性 分 布 的 海拔 最 高 的 极限 环境 ( 吴征 
ARIA EE, 1985; Mani, 1980) 。 这 里 年 均 气 
温 、 积 温 极 低 ， 昼 夜 温差 极 大 ， 强 风 和 日 照 辐 射 
强烈 ， 由 于 冰川 的 剧烈 作用 及 寒 冻 强烈 风化 ， 地 
表 岩 石 剥 落成 大 小 不 等 的 砾石 ， 石 除 中 有 少量 的 
土壤 ， 基 质 贫 将 等 ， 植 物 生长 周期 极其 短暂 ， 也 
最 为 独特 ， 是 环境 条 件 最 为 恶劣 的 高 山 生态 系 统 
( 李 渤 生 等 ,1981; 刘 伦 辉 等 ,1984) 。 由 于 高 山 
冰 缘 带 紧 靠 高 山 永 久 雪 线 ， 环 境 条 件 恶劣 是 高 等 
植物 生存 的 极限 ， 且 特有 性 极 高 ， 是 一 个 自然 和 
相对 独立 的 特殊 生态 系统 。 此 外 由 于 大 部 分 高 山 
山峰 彼此 间 是 独立 的 似 “ 岛 状 ”， 相 应 的 冰 缘 植 
被 带 也 呈现 出 了 “ 岛 状 ”"， 而 使 得 高 山 冰 缘 带 的 
植物 区 系 大 部 分 处 于 “隔离 ”或 “岛屿 ”分 布 
状态 。 在 这 样 极端 严酷 和 高 选择 压力 的 环境 中 ， 
植物 种 类 所 表现 出 来 的 遗传 和 谱系 地 理 结构 以 及 
第 四 纪 冰 期 影响 ， 无 疑 对 探讨 极端 环境 下 植物 多 
样 性 的 形成 分 化 ， 全 面 了 解 高 原 植物 区 系 的 演化 
有 重要 的 意义 。 

因此 ， 本 研究 选择 中 国 喜 马 拉 雅 高 山 流 石 滩 
生境 分 布 的 展 形 科 (Labiatae) 特有 的 单 型 属 绢 
参 属 (Eriophyton) 为 研究 对 象 。 绢 参 ( Eriophy- 
ton wallichii) 是 典型 的 异 花 授 粉 的 多 年 生 草 本 植 
物 〈 吴 征 镑 和 李 锡 文 , 1977) ， 生 于 海拔 2700 ~ 
5 030 m 的 高 山 流 石 滩 ; 主要 分 布 于 青藏 高 原 - 横 
断 山 以 及 印度 、 尼 泊 尔 等 高 山 流 石 浴 。 研 究 中 ， 
我 们 选取 了 核糖 体 (nrDNA) 内 转录 间隔 区 
(CITS) 的 片段 进行 测序 分 析 。77TS 片段 由 于 进化 
速率 较 快 ,测序 难度 较 小 的 特点 常常 被 用 于 较 低 
分 类 阶 元 的 系统 发 育 关 系 研 究 ， 近 年 来 也 较 多 的 
被 应 用 于 推测 物种 谱系 历史 的 研究 ( Collevatti 
等 , 2009; Koch 等 , 2006; 刘 腾 靓 , 2009; 段 义 忠 
等 , 2010) 。 因 此 本 研究 希望 能 揭示 该 种 植物 单 
倍 型 的 地 理 分 布 式样 、 居 群 遗传 结构 等 ， 进 而 推 
测 其 可 能 的 谱系 历史 、 第 四 纪 冰 期 对 其 的 影响 。 
































1 材料 和 方法 
1.1 实验 材料 
研究 材料 采 自 云南 西北 部 、 四 川西 部 、 青 海 和 西 
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藏 ， 履 盖 了 该 种 主要 的 分 布 区 。 按 居 群 采样 标准 对 绵 参 
所 处 地 理 分 布 范围 进行 取样 ， 即 每 个 居 群 随机 采样 6 ~ 
10 个 个 体 ， 每 个 居 群 内 采样 个 体 之 间 至 少 间隔 10 m, FE 
取得 绵 参 20 个 居 群 187 个 个 体 ( 表 1) 。 野 外 采集 生长 
良好 ， 无 病 斑 的 幼 叶 ， 用 硅胶 快速 干燥 后 带 回 实验 室 。 
凭证 标本 存 于 中 国 科学 院 昆明 植物 研究 所 标本 馆 
(KUN), 

1.2 DNA 提取 和 PCR 扩 增 

采用 改进 的 CTAB 法 ( Doyle 和 Doyle, 1987) 从 硅胶 
干燥 的 叶片 中 提取 总 DNA, 用 ITS4 和 ITS5 通用 引物 
(White 等 , 1990) 对 nrDNA 的 ITS 片段 进行 PCR 扩 增 。 
PCR 反应 体系 在 20 pL 体积 中 进行 ,包括 模板 DNA 10 ~ 
40 ng, 1.8 pL 10xPCR 反应 缓冲 液 ，1.4 uL MgCl, (25 
mmol- L^!) , 1.6 uL dNTP 混合 液 (10 mmol L'! ) , IER 
引物 各 1 pL (10 pmol: L), Taq DNA 聚合 酶 0.5 单位 。 
反应 程序 如 下 : JS 94°C 条 件 下 变性 4 min, BU 
94C 1 min, 55C 1 min, 72C 1 min, 循环 35 次 ， 最 后 
在 72% 条 件 下 延伸 7 min。 扩 增 得 到 的 PCR 产物 在 1% 
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1.3 数据 分 析 

对 187 个 个 体 的 PCR 扩 增 片段 的 全 序列 测序 ， 用 
CLUSTAL X 1.81 (Thompson 等 , 1997) 软件 对 测定 的 序 
列 进 行 对 位 排列 ， 并 加 以 手工 校对 。 用 DnaSP 4.0 (Ro- 
zas ^. 2003) 程序 确定 nrDNA 单 倍 型 。 根 据 每 个 居 群 内 
不 同 单 倍 型 的 个 体 数 目 和 总 的 个 体 数 来 计算 单 倍 型 频 
率 ， 然 后 基于 每 个 居 群 的 单 倍 型 组 成 、 不 同 单 倍 型 在 每 
个 居 群 的 个 体 数 目 和 单 倍 型 序列 变异 位 点 ， 分 别 计算 它 
们 的 遗传 多 样 性 大和 (分 别 为 居 群 内 和 居 群 间 平 均 
遗传 多 样 性 ) 以 及 居 群 遗传 分 化 Gy 和 Nop {El (Raymond 
和 Rousset, 1995), (KHE Pons 和 Petit (1996) 所 描述 的 
Jrik, HI PERMUT 软件 计算 。 居 群 间 的 Nsy 值 用 1000 次 
换 进行 置换 检验 。 

应 用 ARLEQUIN 软件 包 version 3. 01 (Excoffier 等 ， 
2006) 中 的 AMOVA ( Analysis of Molecular Variance) 分 
析 方 法 分 别 检测 材料 在 分 布 区 中 居 群 内 和 居 群 间 的 遗传 
变异 水 平 ， 以 及 单 倍 型 分 布 的 Fy 评价 (Weir 和 Cocker- 
ham, 1984) ， 以 进一步 分 析 居 和 群 的 分 化 程度 (1000 KE 







































































] BIOTEKE 公司 小 量 PCR 





的 琼脂 糖 胶 上 电泳 检测 ， 使 月 
产物 纯化 试剂 盒 依照 说 明 书 进 











行 纯化 ， 最 后 在 ABI 3700 








DNA analysis system 测序 仪 上 测序 。 


表 1 
Table 1 


换 检 验 )。 同 时 , 利 




















] ARLEQUIN 软件 包 中 的 Mantel 统 


计 学 检验 (Mantel, 1967) , 分 别 比较 绵 参 居 群 平均 遗传 





距离 与 地 理 距 离 (经 纬度 ) 之 间 的 相关 愧 








from nrDNA haplotypes (H1-H19) in 20 populations of Eirophyton wallichii 





BS 20 个 居 群 材料 来 源 及 nDNA ITS 单 倍 型 (HIL-H19) 的 遗传 多 样 性 (Ad) 组 成 和 频率 


Origin of materials, measurements of haplotypes diversity (Hd) and composition/frequency 
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E (r 值 表示 这 




















































































































编号 (ue. oo e sm mu e 
No. 采集 地 Locality a No. Latitude Longitude Altitude Haplotype 样 指数 样 性 指数 
(N) (E) /m composition Hd T 
1 ”云南 中 多 Zhongdian, YN ZhD 10 28°35’ 99°50’ | 4450 HI15(10) 0 0 
2 WWE Deqin, YN DQ 10 28°20’ 99°04’ 4604 H4(10) 0 0 
3 PYJI| $ 3X Xiangcheng, SC XC 10 29°06’ 100°01' 4790 H16(4); H17(6) 0.533(0.095) 0.00359 
4 四川 理 塘 Litang, SC LT 10 30°10’ 99°59’ 4 420 H10(10) 0 0 
5 ”四 川 牙 江 Yajiang, SC YJ 10 30°00’ | 100*51' 4405 HI2(10) 0 0 
6 “四川 九龙 Jiulong，SC JiuL 10 29°19’ | 101?29' 4581 H11(4); HI2(6) 0.533(0.095) 0.00359 
7 j 藏 左 贡 Zuogong, XZ ZG 10 29°43’ 98°03’ 5028 H18(10) 0 0 
8 i 藏 昌都 Changdu, XZ CD 8 31°18’ 97°32’ 4820 H6(8) 0 0 
9 i 藏 类 乌 齐 Leiwuqi, XZ LWQ 10 31°19’ | 96°30’ | 3965 H7(2); H8(8) | 0.356(0.16) 0.0006 
10 [ 藏 加 查 Jiacha, XZ JC 9 29°02’  Á92?72]!C 4832 H14(9) 0 0 
11 6 藏 错 那 Cuona, XZ CN 7 27°58' 92°01’ 4478 H9(7) 0 0 
12 ”西藏 当 雄 Dangxiong, XZ DX 8 30°36’ 91°06’ 4856 H2(8) 0 0 
13 6 藏 南 木 林 Nanmulin, XZ NML 10 29°15" 87°11’ | 4370 H13(10) 0 0 
14 PREMA Nielamu, XZ NLM 9 28°02’ | 85°57’ 4199 H19(9) 0 0 
15 ”西藏 吉隆 Jilong, XZ JL 10 28933! 85°17’ 4275 H1(10) 0 0 
16 ”青海 称 多 Chengduo, QH ChD 9 33°19’ 97°14’ 4509 H19(9) 0 0 
17 “青海 玉树 Yushu, QH YS 8 32°45’ | 97*12' 4315 H19(8) 0 0 
18 “青海 杂 多 Zaduo, QH ZaD 9 32°57’ 97°11’ 4760 H19(9) 0 0 
19 “青海 曲 麻 莱 Qumalai, QH QML 10 34°09’ 95°56’ 4820 H19(10) 0 0 
20 “四川 得 荣 Derong, SC DR 10 28°48’ 99°12’ 4300 H3(8); H5(2) 0.356(0.16) 0.0006 
合计 187 0.908(0.01) 0.016 
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两 个 矩阵 间 的 相关 性 系数 ) , 并 进行 10 000 次 重复 的 显 
车 性 检验 (Smouse 等 ,1986 ) 。 以 近 缘 种 假 糙 苏 ( Para- 
phlomis javanica) 作为 外 类 群 (Scheen 等 , 2010) ， 应 用 
Network 4.2.0.1 (Forster 等 , 2004) 对 绵 参 20 个 居 群 的 
ITS 进行 基于 简约 性 原则 的 中 间 连 接 网 法 分 析 ( Median- 
joining networks) ( Bandelt &&, 1999), ， 构 建 单 倍 型 网 状 
图 。 同 时 ， 采 用 PAUP* 4.0b10 (Swofford, 2003) 软件 以 
最 大 简约 性 分 析 法 (maximum parsimony, MP) 分 析 各 单 
倍 型 的 系统 发 育 关系 。 分 析 时 所 有 空位 (gap) 作为 缺失 
处 理 , 采用 启发 式 搜索 ，1 000 次 随机 加 入 序列 ，TBR 
枝 长 交换 ， 每 步 保 留 10 棵 树 。 得 到 的 一 致 性 系统 树 分 支 
的 可 靠 性 使 用 “ 靳 带 ”( Bootstrap) 分 析 ， 用 1 000 次 
复 来 检验 单 倍 型 各 分 支 的 支持 率 值 。 

物种 水 平 上 的 居 群 间 的 平均 基因 流 (Nm) 、 单 倍 型 多 
样 性 指数 (Hd) 、 核 苷 酸 多 样 性 指数 (m) (Nei, 1987) 以 及 
Tajima's D (Tajima, 1989) 以 及 Fu 和 Li's D* (Fu FI Li, 
1993) 两 种 无 限 突变 位 点 模型 的 中 性 检验 方法 检测 及 失 配 
分 析 ， 又 称 歧 点 分 布 (Mismatch distribution) 分 析 ( Rog- 
ers 和 Harpending, 1992) 都 在 DnaSP 4.0 程序 中 完成 。 






























































2 结果 
2.1 序列 变异 
对 20 个 居 群 187 个 个 体 进行 了 核糖 体 DNA 
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ITS 片段 的 全 序列 测定 ， 通 过 排序 发 现 ，17S 片 
段 长 度 为 596 bp， 共 检测 到 19 个 单 倍 型 (HI1- 
H19) ( 表 2)， 所 有 单 倍 型 序列 都 提交 到 Gen- 
Bank ， 序 列 号 为 (JN176157-JN176176 ) 。 
2.2 多 样 性 指数 、 和 遗传 结构 分 析 及 失 配 分 析 

在 绵 参 的 20 个 居 群 总 的 单 倍 型 多 样 性 指数 
Hd 为 0.908 (0.01), ER Z FE TETE BL T 为 
0.016。 其 中 单 倍 型 H19 的 频率 最 高 ， 有 45 个 
个 体 拥有 此 单 倍 型 。 每 个 居 群 的 遗传 多 样 性 
(Hd), 、 单 倍 型 数量 、 单 倍 型 组 成 及 频率 见 表 1, 
单 倍 型 的 地 理 分 布 情况 见 图 1 。 

通过 PERMUT 程序 计算 得 出 绵 参 居 群 内 平均 
遗传 多 样 性 五 值 、 总 的 遗传 多 样 性 用, 值 、 居 群 间 
遗传 分 化 Gy 和 Nop 值 分 别 为 : 0.134 (0.0464)、 
0.979 (0.0199) , 0.863 (0.0466) 和 0.957 (0.0171) 。 
利用 U- 统 计 方法 对 绵 参 整个 分 布 区 单 倍 型 变异 的 
地 理 结构 进行 检验 后 发 现 Ny 大 于 Gy 但 不 显著 
(U=2.77) ， 表 明 在 绵 参 分 布 范 围 内 ， 居 群 间 的 遗 
传 分 化 水 平 很 高 (Gop =0. 863) ,但 各 居 群 间 并 无 
明显 的 谱系 地 理学 结构 。 分 子 变 异 分 析 ( AMO- 
VA) 的 结果 表明 ， 居 群 间 的 遗传 变异 为 89.54% , 





























> 人 
Y O i 
e e 
Se 
© 
e, 与 @ 


六 


0 120 240 480 


720 960 








图 1 585 20 个 居 群 中 检测 出 19 个 ITS 单 倍 型 的 地 理 分 布 图 
Fig.1 The geographic distribution of ITS haplotypes (H1-H19) detected in 20 populations of Eriophyton wallichii. 


The pie chart shows the proportions of haplotypes within each population 
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m2 绵 参 核糖 体 I7S 片段 19 个 单 倍 型 间 的 序列 变异 位 点 , 圆 点 (. ) 表示 与 第 一 行 单 倍 型 相同 ,“- ”表示 排序 空位 
Table 2 Variable sites of the aligned sequences of nuclear ribosomal DNA fragments in 19 haplotypes of E. wallichii. 
Dot (. ) indicates that the character states are the same as Haplotype 1, ‘—-’ indicates alignment gap 
1 1 1 1 1 1 1 1 1 22 33 44 44 4 
us 1113444556 6 7 3 4 4 5668990599 0 0 2 2 3 
Haplotype 4 5 9 901356125 0 345626723 0604 3 8 8 9 0 
Hap_1 T CC GA AG AC AG A G GGGGGGGGC-CAG TA TGC NG 
Hap2 A A G G G C AG A A TGT C G 
Hap 3 A A G G G C AG A - C G 
Hap 4 A A G G G C AG A - T C G 
Hap 5 A A G G G C AG A = C G 
Hap_6 - G C C AG A CG C G 
Hap 7 A A G G G C AG C= A CGC = 5s = 
Hap 8 <A A G G G C AG C - ACG 
Hap_9 T GCA C AG A A C G T C G 
Hap_10 G C T C AG A C G A C GA 
Hap_11 G C C AG T A Cc - C G 
Hap_12 G C C AG A C = C G 
Hap_13 G C C AG ACG C G 
Hap_14 G C C AG A C G C G 
Hap_15 G C C AG C G C G 
Hap_16 G C C AG C G C G 
Hap_17 G C C AG C G Gaus = om 
Hap_18 G C AG C - C C G 
Hap_19 G C AG : C = C G 
单 倍 型 4-4 44 4 45 S 5 SS S&S 5S SS 5 SSS FS S S S SH SS F SS SS SS 
33359 900000001 1 1 12223 3 33 4 4 4 4 5 5 5 6 7 
Haplotype 1 2 83353456789%235712315680 37262923 7 
Hapl AGT CA GGCACAAGCGGGGTGTGCGAGGCCTGCGCG 
Hap 2 G G A G G A A C A A G A 
Hap_3 G GA----- - A C A A G 
Hap_4 G GA----- - C A A G 
Hap 5 G GA. G. G C A A G 
Hap 6 TG AG A. G G A C A G A 
Hap7 - - G G , G A AC A G 
Hap_8 G G G A C A G 
Hap_9 G G ; G C T G 
Hap_10 G G ; G A C G A 
Hap 11 G G A G C A A A 
Hap_12 G G A G C A A T G A A 
Hap_13 G G ; G C C A A T 
Hap. 14 G G 2 G C A G A 
Hap. 15 G G : G C A G A 
Hap 16 G G ; G C G A AT 
Hapl7 - - G G , G C G A AT 
Hap. 18 G G j G C T G 
Hap_19 G G . G C T G 


而 居 群 内 的 遗传 变异 为 10. 46% , For = 0. 895 
( 表 3)， 进 一 步 揭 示 了 绵 参 的 遗传 变异 主要 存在 
于 居 群 间 ， 而且 具 有 较 高 的 居 群 分 化 水 平 。 此 











yh, Mantel test 统计 学 检验 比较 了 20 个 居 群 
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示 相 关 性 系数 r=0.18,， P=0.06。 说 明 两 
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者 有 一 定 的 相关 性 ,但 相关 性 很 低 。 说 明 绵 参 各 

居 群 间 没 有 很 明显 的 表现 出 Isolation by Distance 

(IBD) 的 扩散 特征 。 
X3 ”所 有 单 倍 型 的 AMOVA 分 析 

Table 3 Analyses of molecular variance (AMOVA) for all haplotypes 


变异 比例 


























变异 来 源 自由 度 ”方差 变异 成 分 
Anns Variation 
Source of variation df SS VC 
/% 
ie HF Tal . 19 78.279 | 0.42456 Va 89.54 
Among population 
REW . 172 8.533 0.4961 Vb 10.46 
Within population 
总 计 Total 191 86.812 0.47417 
固定 指数 
0.895 


Fixation index F gp 


通过 DnaSP 4.0 程序 在 种 的 水 平 上 对 核糖 体 
DNA 进行 失 配 分 析 ， 得 到 的 失 配 曲线 为 多 峰 曲 
线 (图 2) ， 背 离 了 对 快速 扩张 模型 的 假设 ， 这 
个 结果 在 两 种 基于 无 限 位 点 非 重组 模型 的 中 性 检 
验 (Tajma's D=-0. 19425 (P>0.10), Fu and 
Li's D=2.39723 (P>0.05) 中 得 到 了 支持 ， 表 
明 绵 参 居 和 群 没有 经 历 过 近期 的 扩张 。 
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图 2 绵 参 所 有 个 体 核糖 体 DNA 片段 数据 的 失 配 分 析 
Exp; 期 望 值 ， 0bs: 观测 值 
Fig.2 Mismatch distribution for sequences data of nrDNA fragment 
from all individuals of E. wallichii. The thin represents the 


expectation; The dotted line represents the observed value 


2.3 各 单 倍 型 之 间 的 系统 关系 

用 Network 分 析 得 到 了 各 单 倍 型 网 状 关 系 
(图 3) ， 发 现 除 少 数 几 个 单 倍 型 以 较 少 的 突变 步 
数 直 接 联系 外 ， 大 多 数 单 倍 型 相隔 较 远 。 通 过 众 
多 的 “ 隐 单 倍 型 (可 能 为 未 检测 出 的 单 倍 型 或 
是 已 绝 灭 的 古老 单 倍 型 ) 相连 ， 表 明 现 存 的 各 个 





Paraphlomis 
javanica 





绵 参 19 个 ITS 单 倍 型 的 中 间 连 接 网 状 图 ， 各 单 倍 型 的 饼 图 大 小 对 应 单 倍 型 频率 ， 所 有 突变 步 数 均 已 标 出 ， 
色 圆 点 表示 可 能 未 检测 出 的 单 倍 型 或 是 已 绝 灭 的 古老 单 倍 型 











Fig.3 Median-joining network of 19 ITS haplotypes of E. wallichii. Circumference size is proportional to the haplotype frequency. 


All mutations are marked in the network, black dots represents possible unsampled haplotypes or extinct ancestral haplotypes 
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单 倍 型 具有 相对 较 长 的 独立 发 育 历史 。 而 各 个 
“ 隐 单 倍 型 ”之 间 呈 现 复杂 的 “loop” 网 状 结构 ， 
提示 该 种 植物 历史 上 可 能 有 非常 频繁 的 基因 重组 
事件 。PAUP 最 大 简约 性 (MP) 分 析 ， 得 到 严格 
一 致 性 树 (图 4)。 其 中 聚 在 一 起 且 得 到 较 高 
“Bootstrap ”支持 率 的 单 倍 型 关系 与 Network 分 
析 得 到 的 结果 一 致 ， 如 在 系统 树 上 聚 为 一 枝 的 出 
自得 荣 居 群 的 单 倍 型 H3、H5 以 及 出 自 德 钦 的 
H4 和 当 雄 的 H2， 在 网 状 图 中 也 以 最 短 的 变异 步 
数 连 在 一 起 ， 此 外 来 自 同一 居 群 的 两 种 单 倍 型 如 
类 乌 齐 的 H7, H8 以 及 乡 城 的 H16、H17， 在 系 
统 树 上 和 网 状 图 上 都 呈现 出 最 为 相近 的 关系 。 而 
在 系统 树 上 聚 为 一 枝 的 九龙 居 群 单 倍 型 H11 以 
及 九龙 和 和 雅 江 共享 单 倍 型 H12， 在 网 状 图 中 没有 
直接 相连 ， 各 以 两 步 突变 与 一 个 “ 隐 单 倍 型 ， 
连接 ， 表 明 二 者 具有 各 自 独 立 的 进化 历史 。 
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Paraphlomis javanica 


图 4 绵 参 19 种 单 倍 型 (HI-H19) 的 1 棵 最 大 简约 (MP) 树 ， 
分 支 上 的 数字 表示 简约 性 分 析 1 000 次 重复 得 到 的 靴 带 分 析 
支持 率 (BS) ， 单 倍 型 的 地 理 分 布 情况 标 在 各 单 倍 型 后 


Fig.4 The most parsimonious (MP) tree based on nineteen nrDNA 





haplotypes (HL-H19 ) of E. wallichii, numbers above the branches 
indicate MP bootstrap values (BS) based on 1 000 replicates, and 


geographical distribution of each haplotype was marked out respectively 


3 Wit 
3.1 绵 参 核糖 体 DNA 遗传 多 样 性 和 居 群 结构 

上 述 分 析 表 明 绵 参 的 核糖 体 ITS 片段 总 的 遗 
传 多 样 性 水 平 较 高 ， 有 ,=0.979 (0.0199) 。 远 较 
叶绿体 片段 标记 的 青藏 高 原 的 其 它 植物 高 ， 诸 如 
ROUSE TENE (0.558) (Chen 等, 2008); 长 花 
OG (0.77) (Yang 等 , 2008); eak (0.739) 
(Wang 等 , 2009a) ; 425% MF (0.82) (Wang 等 ， 
2009b) ; 以 及 裸子 植物 祁 连 圆 柏 (0.568) (Zhang 
等 , 2005); 云南 红豆 杉 (0.749) (Gao 等 , 2007) 
和 青海 云 杉 (0.254) (Meng %, 2007) , FA, 
其 居 群 内 的 平均 遗传 多 样 性 Hs 则 显著 小 于 上 述 
几 种 植物 ， 到 =0.134 (0.0464) , fI& Hs 与 高 H, 
表明 绵 参 核糖 体 基 因 组 在 居 群 之 间 分 化 明显 ， 且 
在 绝 大 多 数 居 群 中 单 倍 型 较 单一 。 在 20 个 居 群 
中 ， 其 中 11 个 居 群 每 个 居 群 只 含有 一 种 单独 的 
单 倍 型 有 5 个 居 群 共享 一 个 单 倍 型 ， 此 外 4 个 
居 群 的 单 倍 型 表现 出 多 态 性 ( 表 1)。 

居 群 间 的 强烈 遗传 分 化 也 得 到 来 自 AMOVA 
分 析 的 支持 。AMOVA 分 析 表 明 : 2825 20 个 居 群 
的 mrDNA 遗传 变异 中 ， 居 和 群 间 变异 为 89.54% ， 而 居 
群 内 变异 仅 为 10.46% (363) 。 虽 然 绵 参 核糖 体 
DNA 在 居 群 间 的 分 化 明显 ， 但 是 PERMUT 分 析 
表明 各 居 群 间 并 无 明显 的 谱系 地 理学 结构 。 

高 山 流 石 滩 位 于 自然 树 线 与 永久 雪线 之 间 ， 对 
于 绝 大 多 数 生物 来 说 ， 高 山 流 石 滩 是 陆地 生态 系统 
中 生物 分 布 的 上 限 。 绵 参 的 所 有 生境 仅 限于 高 山 流 
石 滩 ， 在 这 样 极 端的 特殊 生境 中 ， 绵 参 的 遗传 结构 
与 其 分 布 的 特殊 生境 密切 相关 。 青 藏 高 原 高 山 流 
石 滩 呈 现 出 的 孤立 的 “岛屿 ” 状 生境 ， 造 成 了 
绵 参 居 群 间 的 地 理 或 生境 隔离 ， 男 外 较 低 的 种 子 
长 距离 扩散 能 力 ， 限 制 了 居 群 间 的 基因 交流 ， 
此 绵 参 居 群 间 存 在 着 高 度 的 遗传 变异 和 分 化 。 
3.2 绵 参 居 群 的 遗传 结构 演化 

人 们 普遍 认为 ， 在 第 四 纪 冰 期 ,气候 波动 对 
生物 区 系 的 分 布 造成 了 巨大 的 影响 ( Dynesius 和 
Jansson，2000) ， 由 此 可 能 导致 居 群 的 分 化 、 迁 
移 或 灭绝 (Taberlet 等 , 1998) ， 甚 至 加 快 物种 的 
进化 速率 ( Hewitt, 2000; Comes 和 Kadereit, 
1998) 。 如 Qiu 等 (2011) 所 述 ， 基 于 现 有 的 青 
藏 高 原 谱系 地 理学 研究 结果 分 析 ， 青 藏 高 原 植 物 
的 谱系 历史 大 致 可 以 分 为 两 种 情况 。 一 是 一 种 植 
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物 类 群 遵循 “冰期 生境 收缩 ， 冰 期 后 扩张 ” 模 
XX (Hewitt, 2000) ， 高 原 边 缘 低 海 拔 地 区 以 及 横 
断 山 成 为 物种 的 避难 所 ， 保 留 了 大 部 分 的 祖先 单 
倍 型 。 因 此 ， 现 存 植物 在 青藏 高 原 面 上 仅 能 检测 
到 少数 几 种 单 倍 型 ; 二 是 这 些 植物 在 第 四 纪 冰 期 
壬 青藏 高 原 台 面 上 仍 有 部 分 居 和 群 能 够 存活 并 保留 
了 祖先 单 倍 型 ， 青 藏 高 原 面 及 腹地 的 局 部 地 区 成 
为 物种 避难 所 ， 因 而 在 高 原 面 上 也 能 检测 出 多 种 
单 倍 型 。 

绵 参 的 单 倍 型 分 布 和 分 化 情况 则 更 接近 第 二 
种 情况 ， 但 又 有 其 特殊 的 一 面 。 从 绵 参 的 ITS 单 
倍 型 分 布 上 看 ， 在 已 确定 的 19 个 单 倍 型 中 ， 除 
了 单 倍 型 H19 和 H12 由 少数 居 群 共享 外 ， 其 余 
的 17 个 单 倍 型 只 分 布 于 一 个 居 群 内 ， 基 本 上 表 
现 为 “一 个 居 群 一 个 单 倍 型 ”的 模式 ， 类 似 于 

个 个 高 山 孤 岛 的 谱系 地 理 分 布 格局 (图 1) 
居 群 之 间 遗 传 距离 和 地 理 距 离 相 互 关 系 的 分 析 ， 
没有 检测 到 有 意义 的 谱系 地 理 结构 ， 说 明 居 群 之 
间 的 基因 交流 和 联系 较 少 ， 基 于 中 间 连 接 法 
( Median-joining) 构建 的 单 倍 型 网 状 图 也 表明 大 
多 数 单 倍 型 都 位 于 分 支 图 末端 ， 都 是 衍生 单 倍 
型 ， 且 各 单 倍 型 之 间 关 系 较 远 。 因 此 ,与 上 述 物 
种 相 比 ， 绵 参 具 有 一 个 相对 独特 的 谱系 历史 。 显 
然 长 期 的 隔离 应 该 是 造成 这 种 格局 的 主要 因素 ， 
这 也 就 暗示 自 高 山 流 石 滩 生境 形成 后 ， 绵 参 似乎 
就 没有 经 历 过 种 群 分 布 格局 的 变迁 或 交流 ， 也 说 
明 第 四 纪 冰 期 对 该 种 分 布 格局 的 影响 可 能 没有 造 
成 大 的 灭绝 或 迁移 。 

Hewitt (2000) 认为 高 山 带 生境 在 冰期 时 因 
植物 分 布 海拔 下 降 而 相互 连贯 ， 因 而 分 布 较 广 ， 
而 在 间 冰 期 温度 回升 时 ， 植 物 分布 海 拔 上 升 造成 
一 些 物种 发 生生 境 片段 化 。 对 于 绵 参 这 样 的 高 海 
拔 分 布 植物 而 言 ， 现 存 植物 的 间 冰 期 生境 与 冰期 
时 相 比 可 能 有 大 幅度 的 收缩 。 因 此 ， 在 第 四 纪 冰 
期 时 ， 与 大 多 数 温 带 植物 一 般 遵循 的 冰期 退缩 与 
间 冰 期 或 冰 后 期 扩张 模式 不 同 ， 绵 参 可 能 经 历 了 
冰期 的 扩展 和 间 冰 期 或 冰 后 期 种 群 的 收缩 。 这 种 
模式 被 称 之 为 “ 替 位 避难 所 模式 ”。 它 由 Kropf 
ZE (2003) 在 开展 欧洲 阿尔 卑 斯 山 同 域 分 布 但 在 
海拔 上 梯度 差异 的 两 个 物种 的 比较 谱系 地 理学 研 
究 时 首先 提出 。 由 于 从 时 间 长 度 来 看 ， 冰 期 持续 
时 间 较 间 冰 期 长 得 多 (Winograd 等 , 1997; Ruth- 







































































erford 和 D’ Hondt, 2000) ， 因 此 对 适应 寒冷 环境 
的 高 山 植 物 而 言 ， 理 应 有 更 多 的 时 间 来 扩展 居 群 
规模 并 扩大 其 地 理 分 布 区 。 与 之 相反 ， 在 第 四 纪 
冰期 低 海拔 植物 有 可 能 更 容易 遭受 长 时 间 的 物种 
生境 片段 化 或 地 理 隔 离 的 影响 。 因 此 ， 高 海拔 植 
物 居 群 间 的 遗传 分 化 程度 应 比 低 海拔 植物 为 低 。 
日 本 学 者 在 研究 分 布 于 日 本 群岛 的 高 山 植物 类 群 
时 ,检测 到 较 低 的 居 群 间 遗 传 分 化 ， 他 们 认为 这 
些 植 物 类 群 遵循 冰 后 期 生境 片段 化 模式 导致 物种 
的 地 理 隔离 ， 即 冰期 后 种 群 分 布 区 缩小 (Fujii 和 
Senni, 2006; Ikeda 等 , 2006; Ikeda 和 Setoguchi, 
2007) 。 

然而 ， 绵 参 居 群 间 表 现 出 了 较 高 水 平 的 遗传 
分 化 ， 表 明 居 群 间 基因 交流 较 少 。 各 单 倍 型 之 间 
的 网 状 关系 分 析 也 表明 各 现存 单 倍 型 似乎 经 历 了 
较 长 的 隔离 历史 。 我 们 推测 由 于 青藏 高 原 地 形 的 
高 度 异 质 性 ， 对 绵 参 等 高 海拔 分 布 的 植物 类 群 而 
言 ， 在 冰期 时 ， 由 高 海拔 生境 向 低 海 拔 迁 移 ， 因 
其 耐寒 的 特性 加 之 特殊 的 流 石 滩 基 质 ， 使 其 不 可 
能 向 下 迁移 到 很 低 海拔 段 与 该 种 其 它 居 群 连接 或 
融合 。 因 此 在 冰期 居 群 扩展 时 ， 各 居 群 之 间 仍 相 
互 保持 一 种 隔离 状态 ， 很 少 有 机 会 发 生 居 群 间 的 
基因 交流 。 而 在 冰 后 期 或 间 冰 期 生境 海拔 升 高 
时 ， 更 是 加 剧 了 各 个 居 群 之 间 的 隔离 。 因 此 ， 尽 
管 第 四 纪 冰 期 气候 的 波动 使 绵 参 居 群 在 高 山上 下 
迁移 变化 ， 但 彼此 之 间 始 终 保持 隔离 状态 。 长 期 
的 地 理 隔离 导致 居 群 间 高 度 分 化 ， 并 且 每 一 个 居 
群 固定 了 一 种 单 倍 型 。 绢 参 的 遗传 变异 主要 存在 
于 居 群 间 ， 只 有 少数 几 个 单 倍 型 分 布 于 两 个 以 上 
的 居 群 中 ， 如 最 常见 的 单 倍 型 H19 WER Xx 
DUK. PKA. ARS AH RIESE; 单 倍 型 H12 
为 九龙 和 和 雅 江 居 和 群 所 共享 。 这 些 少数 共享 的 单 倍 
型 可 能 是 由 于 古老 单 倍 型 的 残 遗 ， 也 可 能 由 于 侦 
然 的 长 距离 扩散 或 基因 交流 形成 。 

失 配 分 析 结 果 为 多 峰 曲 线 ， 表 明 绵 参 并 未 在 
冰期 后 经 历 过 种 群 的 快速 扩张 。 我 们 的 研究 结 
与 先前 研究 过 的 青藏 高 原 植物 不 同 ， 绵 参 呈现 出 
高 度 分 化 的 遗传 结构 ， 常 常 呈现 出 “一 个 居 群 ， 
一 种 单 倍 型 ”的 模式 ,这 可 能 是 由 于 第 四 纪 冰 
期 气候 波动 以 及 青藏 高 原 极端 地 形 的 组 合 效 应 。 
是 否 还 有 其 它 流 石 滩 植物 表现 出 相同 的 遗传 结构 
有 待 进一步 研究 。 
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